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An extension of the Ambegaokar-Baratoff relation to a superconductor-insulator-superconductor 
(SIS) Josephson junction with multiple tunneling channels is derived. Appling the resultant relation 
to a SIS Josephson junction formed by an iron-based (five-band) and a single-band Bardeen-Cooper- 
Schrieffer (BCS) type superconductors, a theoretical bound of the Josephson critical current (/c) 
multiplied by the resistance of the junction {Rn) is given. We reveal that such a bound is useful 
for identifying the pairing symmetry of iron-pnictide superconductors. One finds that if a measured 
value of IcRn is smaller than the bound then the symmetry is zbs-wave, and otherwise s-wave without 
any sign changes. In addition, we stress that temperature dependence of IcRn is sensitive to the 
difference of the gap functions from the BCS type gap formula in the above heteroj unction. 

PACS numbers: 74.50.+r,74.70.Xa 



I. INTRODUCTION 

Since the discovery of iron-based superconductors^^, 
their pairing symmetry has been intensively debated. 
According to the spin fluctuation mechanism associated 
with the Fermi surface nesting, ±5-wave symmetry was 
proposed as a pairing scenario However, the debate 
has not been settled down. The ±5- wave symmetry is ex- 
pected to be fragile against non-magnetic impurities 
Some experiment supported this idea, while the oth- 
ers—""— presented controversial results. Hence, a direct 
and unambiguous evidence like the phase sensitive mea- 
surement in High-Tc cuprate superconductors^^ is now 
in great demand. In fact, a large number of the methods 
to seek a definite signature have been examined, e.g., 
tunneling spectroscopy—"—, corner junctions——, ob- 
servation of half-integer flux-quantum jump^^, scanning 
tunnel microscopy and so on. 

Josephson junctions are sensitive devices reflecting su- 
perconducting states of each electrode. Very recently, 
various types of Josephson junctions with iron-pnictide 
superconductors were successfully fabricated and typi- 
cal Josephson effects were conflrmed— ~— . Among them, 
a Josephson junction between an iron-based and a con- 
ventional 5- wave single-gap superconductors has been re- 
garded as a possible candidate to directly detect the pair- 
ing symmetry of iron-based superconductors. The het- 
erojunction system is theoretically described by multiple 
tunneling channels, some of which are tt channels and 
the others are ones— i^S. Authors suggested anomalous 
critical current reduction Riedel anomaly cancella- 
tion^^, and enlargement of the Josephson vortex core^^. 

In this paper, we derive Ambagaokar-Baratoff rela- 
tion— in the heterojunctions with multiple tunneling 



channels and clarify that a theoretical bound of IcRn 
products distinguishes ±5- wave from 5- wave without any 
sign changes, which is simply denoted as 5- wave through- 
out this paper. We examine two kinds of materials, 
(Ba,K)Fe2As2 (122 compound) and LaFeAs(0,F) (1111 
compound) as the iron-based superconducting electrode 
in the heteroj unction. Employing the density of states 
(DOS) ratios and the superconduting gap ratios given by 
flve-band quasi-classical theory with the flrst-principles 
calculation^^, the theoretical bounds are evaluated. The 
temperature dependences of IcRn are also demonstrated 
in both ±5- wave and s-wave. 

The paper is organized as follows. Section [III is the 
derivation of the Ambegaokar-Baratoff relation in the 
junction with multiple tunneling channels. Based on the 
result, we propose a criterion for identifying the pairing 
symmetry of iron-based superconductors. The key crite- 
rion is an upper bound of the Josephson critical current 
for the its-wave, which corresponds to a lower bound for 
the s-wave. In Sec. lIIH we apply this criterion to typi- 
cal iron-pnictide superconducting materials and theoret- 
ically conflrm its effectiveness. Section [IVl is devoted to 
the summary. 



II. THEORETICAL BOUNDS OF JOSEPHSON 
CRITICAL CURRENTS 

We examine a superconductor-insulator-superconductor 
(SIS) Josephson junction, as shown in Fig.[Tl The elec- 
trode 1 (2), whose length is s' (s) in the direction of the 
z axis, is a single-band (flve-band) superconductor. The 
insulator, whose length is d in the direction of the z axis 
and the dielectric constant is e, is sandwiched between 
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FIG. 1: (color online) Schematic diagram of a SIS Josephson 
junction between five-band and single-band superconductors. 
Here the electrode 2 is assumed to be iron-based supercon- 
ductor. The three superconducting phases enclosed by the 
solid line can be assigned as hole bands of iron-based super- 
conductor, while the two phases enclosed by the dashed line 
as electron bands. 



the two different superconducting electrodes. 

One of the fundamental quantities characterizing 
Josephson junctions is the Josephson critical current 
density. One can find various discussion for the cases 
including two-band superconductors in several refer- 



ences- 



25.39,44 46 



. However, there is no work on arbitrary 
A/'-band superconductors. This paper treats such a gen- 
eral case. Throughout this paper, we assume that there 
is no relative superconducting phase fluctuation in multi- 
band superconducting electrodes. In this case, the sys- 
tem is described by an electric circuit as shown in Fig. [21 
We remark that such a rigid parallel circuit modeling is a 
good description as far as the relative phase fluctuations 
are fully pinned. We will discuss corrections of relative 
phase fluctuations to Josephson effects elsewhere. Un- 
der this assumption, the total Josephson critical current 
density is given by 



E(ii) 
Jicosxo • 



(1) 



(a) 

Xo is a constant phase between the ith and the jth 
superconducting gaps, which reflects the symmetry of a 



static gap solution. Equation ([T]) includes tt channels 
when sign changes occur between the superconducting 
gaps. The basic formalism to derive Eq. ([T]) is shown 
in AppendixlAl For the s-wave case, Xo*^^ = for any 
i. On the other hand, a part of {Xo^^} should be tt 
for the ±5- wave symmetry. As an example for the ±s- 
wave case in the electrode 2, we take Xo^^^ = Xo^^^ = 

Xo^^^ = and Xo^^^ = Xo^^'' = tt as schematically shown 
in Fig. [21 Equation ([T]) indicates that jc for the ±5- wave 
symmetry is always smaller than that for the s-wave, i.e., 
jc(s-wave) > jc(±5-wave) 

Let us turn to a microscopic formula for ji in Eq. ([T]) . 
First, we give notations. As for the electrode 1, we denote 
the DOS on the Fermi surface and the superconduct- 
ing gap amplitude as, respectively, A^s and A*^^^(> 0). 
Similarly, as for the electrode 2, Ni and A^*^(> 0) are 
the ith DOS and the ith gap amplitude, respectively. 
In addition, we define "smaller" and "larger" gaps as 
As,i = min{A(*), A(^)} and AL,i = max{A(^), A^^)}, re- 
spectively. Thus, assuming the full gap solutions in both 
superconducting electrodes, we microscopically calculate 
ji using a standard second order perturbation theory 
with respect to a tunneling channel Then, we have 



Ji 



1 1 TrAe 



2e 



where 

1 47re2 



h 



and W is the area of the junction interface. In the defi- 
nition of r~], the tunneling constant associated with the 

ith tunneling channel is denoted as T*^*^ . The quantity ki 
corresponds to the ratio of the smaller gap to the larger 
one, ki = [1- (As,i/AL,i)^]^/^ The function K{k;iy) is 
given by 



K{k;u) 



I 



^ tanh(i/Vl - A:2xV2) 

y^{l-k^X^){l-X^) 
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FIG. 2: Electric circuits corresponding to a SIS Josephson 
junction. The circuit has the five parallel branches associated 
with the tunneling channels. The bias current is denoted as 
^bias- The left panel corresponds to the case of the s-wave, 
while the right panel to the it s-wave. 



Combining the above arguments with Eq. ([T]) , we ob- 
tain 



IcRn 



Rn 7rAeff,i (il) 

^ cosxo , 



2e 



(2) 



where Ic = jcW. The combined resistance = 
experimentally measured when a bias 
current is greater than /c, while the individual measure- 
ment of rn,i is practically impossible. Equation (|2]) is a 
generalized formula of the Ambegaokar-Baratoff relation 
for multi-channel hetero junctions. Brinkman et at. and 
Agterberg et al. obtained similar results in the context 
of MgB2, i.e., two-band superconductor. 

Here, we derive a simple relation from Eq. (j2|) when 
the electrode 2 has the s-wave symmetry. The resultant 
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Ni/Ntot 


iV2/iVtot 


Ns/Ntot 


iV4/iVtot 


m/Ntot 


BaFe2As2 


0.1687 


0.2757 


0.2104 


0.1803 


0.1649 


LaFeAsO 


0.1133 


0.3009 


0.2323 


0.1608 


0.1928 



TABLE I: DOS ratios on the Fermi surfaces of iron-based materials evaluated by a first-principles calculation—. The first three 
DOS's (A^i, A^2, and N3) corresponds to the hole bands, while the remaining DOS's (A^4 and N5) to the electric bands. 



(T = OK) 




A^2^ /A 


A^^^ /A 


A^^^ /A 

^ / ^max 


A*^^^ /A 

^ / ^max 


2An.ax/feTr" 


A*^^^ /A 

^ / ^max 


(Ba,K)Fe2As2 


0.9395 


1 


0.5189 


0.9415 


0.9691 


3.785 


0.192 


LaFeAs(0,F) 


0.5052 


1 


0.2677 


0.5084 


0.5300 


4.426 


0.232 



TABLE IL Superconducting gap amplitude ratios of iron-pnictide materials at zero temperature estimated by a five-band 
quasi- classical theory—. 
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FIG. 3: (color online) Temperature dependence of iron- 
pnictide and single-band superconductor gap amplitudes. As 
for the iron-pnictide materials, we employ the previous re- 
sults The gap amplitudes are normalized by eVb = 
e[A(^)(0) + A™ax(0)]. The solid lines (A^^^ A^^), and A^^^) 
correspond to the hole band, while the dashed lines (A*^^-* and 
A^^-*) the electric band. In the case of the =bs-wave symme- 
try, the gaps corresponding to the dashed lines have the rel- 



ative minus signs (i.e., Xo 



(41) 



(Te^ 



(51) 

Xo 



: 38 K) and (b) LaFeAs(0,F) {T^ 



tt). (a) (Ba,K)Fe2As2 



27 K). 



expression can give useful information about the pairing 
symmetry. The right hand side of Eq. (j2j) is a summation 

of positive quantities in this case (i.e., Xo*^^ = foi" 
i). Accordingly, we find that 



/c(5-wave)i?n ^ 



ttA^ 



A* = min Aeff,i, 



(3) 



which gives a theoretical lower bound of IcRn for the 
5- wave. From the above argument, we find that if the 
symmetry of the electrode 2 is 5- wave then Eq. (|3]) must 
be fulfilled. Namely, if a measured value of IcRn satisfies 



the inequality 



IcRn < 



ttA* 
2e 



(4) 



then one can exclude the possibility of the 5- wave. This 
method does not give more detailed information about 
the symmetry but clarify whether tt channels exist in 
the multiple tunneling junctions with iron-pnictides. The 
present scheme, "lower-bound criterion" is a simple and 
convenient classification of iron-based superconducting 
materials. 

Let us discuss experimental applicability of the crite- 
rion. The direct experimental data for bulk supercon- 
duting samples, e.g., angle resolved photoemission spec- 
troscopy— 1^, provides the gap amplitudes in the 1-gap 
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FIG. 4: (color online) IcRn products normalized by Eq. (|5}. 
The red squares are for the s-wave, while the green crosses 
are for the =bs-wave symmetry (i.e., Xo^^'* — Xo^^ — ^ind 
Xo^^'* = Xo^^'* — We also show the theoretical lower bounds 
for the s-wave symmetry at zero temperature (7rA*(0)/2eVb), 
depicted as the blue solid lines. We find that IcRn/Vo for zbs- 
wave is much smaller than the theoretical lower bound for the 
s-wave at zero temperature, (a) (Ba, K)Fe2As2 (T^ = 38 K) 
and (b) LaFeAs(0,F) {T^^^^ = 27 K). 
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FIG. 5: (color online) Results for smaller iron-pnictide superconducting gap amplitudes. The DOS and the gap ratios are 
the same as the ones in Fig.SJa), but the iron-pnictide superconducting transition temperatures are much smaller than in the 
previous example. The IcRn products normalized by Eq. ([5]) are shown in (c) ((d)), in which the gap amplitudes drawn in (a) 
((b)) are employed. 



and the iron-based superconducting electrodes. There- 
fore, in principle, one can input A* in Eq. (|3]). How- 
ever, super conduting gap amplitudes relevant to Joseph- 
son junctions is generally smaller than the ones in the 
bulk superconductors due to the damage piled up in in- 
terface fabrication. It indicates that direct comparison of 
a measured IcRn product to A* based on the bulk data 
may give no practical information. Therefore, we focus 
on a transition voltage of the all tunneling channels into 
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FIG. 6: (color online) Comparison of iron-pnictide's gap func- 
tion to the single-band BCS gap function. The cyan square is 
the small gap A*^^-* calculated by the five-band quasi- classical 
theory for (Ba,K)Fe2As2 while the black solid line corre- 
sponds to Eq. (|6|), in which we use T^^^ instead of T^. One 
finds that a difference between the two functions exists around 
^^^^122 ^ 0.4-0.8. We can find a similar discrepancy for A^^^ 
which is much smaller than in the present figure. 



^0 = -[A(^ni^ = 0) + A^ax(T = 0)], 



(5) 



in which A^ax = max^ A^*^. One can measure the sup- 
pressed Vq experimentally if the junction shows hysteretic 
I-V characteristics^^. We emphasize that a scaled quan- 
tity A*/e]^ is given by only the gap amplitude ratios 
since the bulk gap ratios are expected to be kept as long 
as the damage is not too severe. Therefore, the scaled 
lower bound A*/el^ is experimentally evaluated by the 
gap amplitude ratios for the bulk samples, and this quan- 
tity should be compared to measured IcRn/VQ- Conse- 
quently, we conclude that the inequality (H]) is an effective 
formula to examine the pairing symmetry. 



III. APPLICATION OF LOWER-BOUND 
CRITERION TO IRON-PNICTIDE 
SUPERCONDUCTORS 

Let us evaluate the right hand side of Eq. ([2|) in real 
iron-pnictide materials and check how the criterion based 
on the inequality (j4]) works. For the sake of simplicity, 
we employ a simple model for the tunneling constants. 
Namely, we assume that T^*^'s take channel-independent 
constants. It means that Rn/rn,i is equal to A^^/A^tot, 
where A^tot = Xli ^i- Then, we find that IcRn/Vo is sim- 
ply a function of the superconducting gap ratios (i.e., 
A^'V^max and A^^V^max) and the DOS ratios (i.e.. 

We concentrate on two iron-pnictide superconduct- 
ing materials, (Ba,K)Fe2As2 (T^^^^ = 38 K) and 
LaFeAs(0,F) {T^^^^ = 27K). Here, we denote each su- 
perconducting transition temperature as T^'^'^ or T^^^^. 
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As for the single-gap superconducting electrode mate- 
rial, we choose an alloy of Pb, i.e., Pb-In-Au. The 
superconducting transition temperature = 7K and 
A(^)(0)//cbT^' = 1.98, respectively^^. The gap am- 
plitude ratios and their temperature dependence of 
(Ba, K)Fe2As2 and LaFeAs(0, F) were evaluated by five- 
band quasi-classical theory combined with the DOS's via 
the first-principles calculations and several kinds of ex- 
perimental data for the bulk properties We label three 
hole bands as 1 to 3 and two electron ones as 4 to 5. The 
sign changes are assumed to occur between the hole and 
the electron bands. The DOS and the gap ratios are sum- 
marized in Tables |T] and |TT1 respectively. As for the tem- 
perature dependence of the single-band superconducting 
gap, we utilize the BCS type gap formula — 



A(^)(T) = A(^)(0)tanh<^ A 



B{ ^ 

T 



(6) 



where A = 1.82, B = 1.018, and C = 0.51. 

Figure [3] displays the temperature dependence of the 
superconducting gap amplitudes. Figure |4] shows the 
temperature dependence of IcRn normalized by Vq. The 
lower bound at T = given by the right hand side of the 
inequality (|3]) is also depicted by the (blue) horizontal 
line. Due to jc(s-wave) > jc(±5-wave), IcRn for the ±5- 
wave becomes smaller than the one for the s-wave over 
all temperature regions. Here, let us focus on the zero 
temperature. IcRn/Vo for the ±s-wave is smaller than 
the lower bound (7rA*(0)/2eVb) while that for 5-wave is 
larger for both of the iron-pnictide superconducting ma- 
terials. Hence, the comparison of a measured value of 
IcRn to the theoretical bound gives a useful criterion for 
the symmetry in the iron-based superconductors. 

Next, we investigate specific cases in which a part of 
or all iron-based superconducting gaps are smaller than 
the single-band BCS gap, as shown in Fig.[5fa) and[5fb). 
In such a study, we keep the gap ratios to be the val- 
ues for (Ba,K)Fe2As2 shown in Tables |T] and |TT1 except 
for A*^®^/Amax- Figures [5fc) and [5fd) show the IcRnS 
for T^^^^ = 10 K and 6 K, respectively. We find again 
that the inequality (j4]) is fulfilled at zero temperature. 
It means that the present criterion works in every case. 
On the other hand, we notice that the temperature de- 
pendences of IcRn are relatively anomalous. We find a 
convex behavior in the middle temperature range (i.e., 
2^^ja22 ^ 0.2-0.8) contrary to our naive expectation. 
Such peculiarity comes from a discrepancy between the 
gap functions calculated by the five-band quasi-classical 
theory and the BCS type gap formula (j6|), as shown in 
Fig.O In the previous case as Fig. HI the single-band BCS 
gap solely provides the contributions to the temperature 
dependence of IcRn- However, the present case refiects 
all the gap temperature dependences. Thus, the temper- 
ature dependence of IcRn is found to be quite sensitive to 
the difference of the gap functions from the single-band 
BCS type gap formula—. 



IV. SUMMARY 

We derived the Ambegaokar-Baratoff relation in the 
SIS Josephson junction with multiple (more than two) 
tunneling channels and proposed a criterion to identify 
the pairing symmetry of the iron-based superconductors. 
If a measured value of IcRn product is smaller than the 
lower bound for the 5- wave, then one concludes that the 
symmetry of the superconducting electrode is the ±5- 
wave. We actually revealed that the criterion well works 
in the typical iron-pnictide superconductors by employ- 
ing the DOS and the gap ratios calculated by the five- 
band quasi-classical theory and the first-principles calcu- 
lation. In addition, the theory predicted that the temper- 
ature dependence of IcRn is sensitive to the deviation of 
the temperature dependence of the gap from the single- 
band BCS formula. The method is simple and conve- 
nient in contrast to the phase sensitive measurement like 
TT junctions which require much more elaborate setups. 

The Ambegaokar-Baratoff relation is one of the fun- 
damental identity in the SIS Josephson junction. The 
present analysis suggested that this type of the relation 
provides more fruitful information when the junction has 
multiple tunneling channels. 
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Appendix A: Formalism for gauge- invariant phase 
differences 

A basic formalism for the hetero Josephson junction 
dealt in this paper is presented. Assuming uniformity 
along y axis, the effective Lagrangian density on the zx 
plane— is given by 



87r/i' 



■^^cosx(^'^) + |^(^ll)^ (Al) 



i<i' 



Stt 



^{i) _ (^(s)(e*d/nc)A|i and x^''^ = ^^'^ 



(^vv = 6>(* ) — 0^^\ X^* is the relative phase differ 



where 0^'^^ 

ence between the different superconducting gaps. The 
first and the second terms in Eq. ()Aip represent charge 
compressibility in the electrode 1 and 2, respectively. 
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where = (n/e*)a^(/?^^) + A?, a? = {h/e'')dt^^^ + A^, 
and e* = 2e. The electric scalar potential in the elec- 
trode i{= 1,2) is denoted as A^, and the charge screen- 
ing length in the electrode 1 (2) is written as ji' {fii). 
The third term in Eq. (|A1|) is the Josephson coupling, 
in which hji/e* is the coupling constant associated with 
the ith tunneling channel. The forth term in Eq. ()Aip 
is called the inter-band Josephson coupling energy 
whose origin is the inter-band interaction between dif- 
ferent bands in the electrode 2. In the gauge- invariant 
phase difference 0^^\ the vector potential in the insula- 
tor = J^^^^^A^{z)dz^ in which A^{z) is the z- 
component of the vector potential. The electric field is 
expressed as E^-^ = -c'^^dtA^-^ - - A^). 

As the Euler-Lagrangian equation with respect to A'^, 
we obtain the Josephson relation 

Y^^d,e(^) = '^E^„ (A2) 

i 

where a' = eji''^ / s' ai = e/if/sd, = J2i^7^ 

A = 1 + a' + q;. Next, as the Euler-Lagrangian equation 

with respect to A^i, we have the Maxwell equation 

= ?ir^*^2i + E ^- (A3) 

Combining Eq. (|A3|) with Eq. (|A2[) leads to the equation 



If one adds a dissipation and an external bias current 
terms to Eq. ()A4p , one has the resistively and capacitively 
shunted junction model with multi-tunneling channels. 

Let us turn to a formula for the Josephson criti- 
cal current density. Generally, the system has a rela- 
tive superconducting phase fiuctuation originating from 
the inter-band Josephson coupling ^^^^f Jn' cosx'''^ 
which generates Josephson-Leggett collective excitation 
modes The equations of motion for x^* can be 
obtained as the Euler-Lagrange equations with respect 
to cp'^^^ and cp'^'^^—. In this paper, we simply assume that 

each x^* is fixed as a constant Xo* refiecting the sym- 
metry of a static gap solution. Such an assumption can 
be validated when \Jii'\ ^ ji,j2- In this case, we have 
0^'\t) = 0^^\t) + Xo^\ where Xo = or tt. Equation 
()A4p is then rewritten by 



^^^tO^''^ + E c°«Xo = Jbias. (A5) 



Here, we add a bias current density jbias to the right hand 
side. Equation ()A5p indicates that a static solution exists 

as long as jbias < I Ei cosxo ^^1- 
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